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ABSTRACT 

There are two fundamentally different physical origins of faint satellite galaxies: cosmo- 
logical sub-structures that contain shining baryons and the fragmentation of gas-rich 
tidal arms thrown out from interacting galaxies during hierarchical structure forma- 
tion. The latter tidal-dwarf galaxies (TDG) may form populations with correlated 
orbital angular momenta about their host galaxies. The existence of TDGs is a strin- 
gent necessity because they arise as a result of fundamental physical principals. Wc 
determine the significance of the apparent disc-like distribution of Milky Way (MW) 
satellite galaxies. The distribution of the MW satellites is found to be inconsistent 
with an isotropic or prolate DM sub-structure distribution at a 99.5 per cent level 
including the recently discovered UMa and CVn dwarf spheroidal galaxies. The dis- 
tribution is extremely oblate and inclined by about 88° with respect to the the MW 
disc. We also apply the methods to Andromeda's (M31) satellite galaxies using two 
recently published data-sets. It can not be excluded that the whole population of 
M31 companions is drawn randomly from an isotropic parent distribution. However, 
two subsamples of Andromeda satellites are identified which have disc-like features. A 
kinematically motivated subsample of eight Andromeda satellites forms a pronounced 
disc-like distribution in both data-sets. The existence of this disc would be inconsis- 
tent with a CDM parent distribution of subhaloes if the disc is rotationally supported. 
The M31 satellite distribution is inclined by about 59° with respect to the M31 disc, 
and has virtually the same orientation as the disc derived for the whole M31 satellite 
sample. We present a new geometric method to set restrictions on possible locations 
of angular momentum vectors for Andromeda satellites. Our conclusion is that both, 
the MW and M31, may indeed have satellite galaxies derived from TDGs. Further, 
both host-discs and both identified discs-of-satellites are highly inclined relative to the 
supergalactic plane. The discs-of-satellitcs therefore cannot be created from individual 
accretion events from the supergalactic plane further supporting the possibility that 
they are of TDG origin. 

Key vi^ords: Galaxies: evolution. Galaxies: formation. Galaxies: structure. Galaxies: 
dwarf, Galaxies: Local Group, Galaxies: fundamental parameters 



1 INTRODUCTION 

It is well known that in cold-dark-matter (CDM) simula- 
tions of large-scale structure in the Universe the growth of 
structure proceeds via a sequence of hierarchical collapses 
driven by gravity. Instabilities, perturbations, and torques 
in the accretion process generate filamentary-like networks 
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that agree well with observed distributions of galaxies from 
large-scale redshift surveys. Systems coalesce onto these fila- 
mentary networks, which in turn merge to form the structure 
we see today. 

There are, however, several issues on smaller scales 
which are difficult to address using CDM simulations. The 
number of observed satellite galaxies around the Milky 
Way (MW) and Andromeda (M31) is significantly smaller 
than the number predicted by models (Moore et al. 1999; 
Klypin et al. 1999; Governato et al. 2004, but see also 
Kase et al. 2006). This so-called 'sub-structure crisis' is usu- 
ally addressed by invoking small-scale baryonic processes or 
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baryonic-dark-matter biases (Kazantzidis et al. 2004, and 
references therein) , but even then the central density profiles 
of the sub-structures remain cuspy, despite tidal heating and 
destruction in the host halo. Even extreme baryon removal 
cannot evolve a cusped to a cored DM halo (Gnedin & Zhao 
2002). The cuspy profile is in disagreement with the density 
profiles inferred for well-observed MW dSph satellites that 
are interpreted to be the most dark-matter dominated ob- 
jects known (Wilkinson et al. 2002; Kleyna et al. 2003), an 
interpretation that may need revision (Munoz et al. 2005). 
Indeed the DM profiles inferred for the dSph satellites by 
solving the Jeans equation are completely inconsistent with 
CDM profiles (Wilkinson et al. 2006; Gilmore et al. 2006). 

The spatial distribution of the Milky Way satellite 
galaxies has long been known to show asymmetric pat- 
terns and probable streams of satellites (e.g., Kunkel & De- 
mers 1976; Lynden-BeU 1976; Majewski 1994; Lynden-BeU 
& Lynden-BeU 1995; Hartwick 2000, Palma, Majewski & 
Johnston 2002). Kroupa, Theis & Boily (2005) tested the 
spatial distribution of the satellite system against the null- 
hypothesis that it is drawn randomly from a spherical dis- 
tribution of dark-matter dominated subhaloes. They found 
that this hypothesis can be excluded with very high statis- 
tical significance, given that empirical constraints show the 
MW potential to be spherical (e.g., Fellhauer et al. 2006). 

In reply to Kroupa et al., Kang et al. (2005) argued 
that if the Milky Way satellites follow the distribution of the 
dark-matter within the MW halo rather than the distribu- 
tion of substructure selected by present-day mass, then the 
observed distribution of the MW satellites is consistent with 
being CDM sub-haloes. They based their argument mainly 
on the apparent rms-height of the observed disc-like distri- 
bution showing that a steeper radial number-density distri- 
bution yields a smaller rms-height. Zentner et al. (2005) used 
a semi-analytic model to identify luminous satellite galax- 
ies in CDM host-haloes. As a second test they tagged the 
most massive dark-matter sub-haloes as luminous satellites 
(Stoehr et al. 2002). They showed that an isotropic distri- 
bution is not the correct null-hypothesis, but that the host 
haloes are mildly triaxial, tending to be more prolate than 
oblate. Based on the relative height of the distribution they 
argued that the MW satellite system is consistent with be- 
ing CDM substructure, albeit with a low probability. Sim- 
ilarly, Libeskind et al. (2005) identified luminous satellites 
using a different semi-analytic model for star formation. Us- 
ing halo merger trees, they found that the distribution of the 
most massive progenitors is consistent with the observed dis- 
tribution of the MW satellites in all their simulations and 
that a spherical parent distribution is not the correct null- 
hypothesis. In contrast to Zentner et al. (2005) they also 
found that the distribution of the most massive sub-haloes at 
present is significantly different from that of the MW satel- 
lites and the most massive progenitors. All of these results 
offer different solutions to the disc-of-satellites problem, but 
these simulations are based on CDM models that do not in- 
clude the dissipative physics of galactic disc-formation. The 
existence of a disc-galaxy and the orientation of the baryonic 
disc relative to the satellite distribution need to be postu- 
lated. 

Similarly to the Milky Way satellite system, the satel- 
lites of Andromeda seem to be anisotropically distributed 
as weU (Grebel, Kolatt & Brandner 1999; Hartwick 2000). 



Koch & Grebel (2006) addressed this issue by performing 
an analysis using a great-circle fitting routine. They found 
a planar-like distribution with low statistical significance. 
However, for a morphologically motivated subsample, in- 
cluding most dSph/dE satellites, they claim a highly sig- 
nificant polar great plane. McConnachie & Irwin (2006b) 
showed that the M31 satellite system is significantly skewed 
in the direction of the MW. They also identified possible 
ghostly streams (Lynden-BeU & Lynden-BeU 1995) of sub- 
samples of satellites based on the intersection of all possi- 
ble kinematical poles (i.e. directions of orbital angular mo- 
menta) . 

Given the problems the CDM hypothesis has in dealing 
with virtually all aspects of the dwarf-satellite problem, it is 
useful to step-back and to consider some issues of fundamen- 
tal physics: dwarf satellite galaxies can have two fundamen- 
tally different origins (Hunter et al. 2000): either they are 
hierarchical building blocks, DM sub-haloes (in CDM cos- 
mology) , containing shining baryons and are not yet merged 
with a host galaxy (e.g., Read et al. 2006), or they are anti- 
hierarchically formed as tidal-dwarf galaxies (TDGs) in tidal 
arms thrown out from interacting gas-rich galaxies. While 
the former dwarfs critically depend on cosmological theory, 
TDGs are a result of well-established fundamental physical 
principles, the conservation of energy and angular momen- 
tum. TDGs must therefore arise in any cosmological theory 
of structure formation. The formation of TDGs is observed 
in the local Universe (e.g., Hunsberger et al. 1996; Kroupa 
1998; Weilbacher et al. 2003; WaUer et al. 2006), the TDG 
candidates having gas masses of up to 10^ Mq and ongoing 
star formation. The efficiency of the production of TDGs 
is expected to be much higher during early cosmological 
epochs due to the large gas content of the progenitor galax- 
ies. Groups of TDGs originating from one encounter have 
correlated orbital angular momenta and may therefore later 
form a disc-of-satellites. 

In the following the spatial distribution of both, the 
MW and M31 satellites, is investigated using the same anal- 
ysis methods, allowing a direct comparison of the properties 
of the satellite distributions. Our ansatz is to test the null- 
hypothesis that the satellites of the MW and M31 are dis- 
tributed in a disc. Exclusion of this hypothesis would only 
imply probable consistency with the theoretical DM sub- 
structure distribution, notwithstanding the failure of the 
CDM hypothesis to account for the inferred profiles and 
number of the putative DM haloes of the satellites (Gilmore 
et al. 2006). A pronounced disc- like distribution would pro- 
vide strong support for causally-connected satellites, subject 
to the condition that their angular momenta are correlated. 

We describe the mathematical methods used to fit 
planes (§2.1) and to analyse the data (§2.2). We analyse the 
spatial distribution of the Milky Way (§3.1) and Andromeda 
(§3.2) satellite system, respectively, using two recent data- 
sets for the latter. A new method based on radial velocity 
measurements is used to set some constraints on possible 
kinematic associations of Andromeda satellite galaxies (§4). 
Finally, in §5 we constrain the shape of the parent distribu- 
tions of the MW and M31, which may be prolate, spherical, 
or oblate DM haloes, and discuss the results in §6. 
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2 TECHNIQUES 
2.1 Plane fitting 

To fit a plane to tlie data, an unweiglited fitting algoritlim, 
icnown as an algebraic least-squares (ALS) estimate method 
(see, e.g. Chojnacki et al. 2000) or eigenvalue analysis, is 
incorporated. It is similar to the algorithm used by others 
(e.g. Dubinski & Carlberg 1991; Hartwick 2000; Libeskind 
et al. 2005) . In this method the centroid of the data points, 
ro, is calculated and an eigenvalue analysis of the moment 
of inertia tensor To of the position vectors Vi = Vi — ro, 
i = 1 . . . n, n being the number of satellites, relative to the 
centroid ro is performed. The eigenvector corresponding to 
the smallest eigenvalue is the normal of the plane and the 
plane contains the centroid. Incorporating the centroid of 
the data ensures that we correctly find the plane that has 
the minimum orthogonal distance to the satellites. Not con- 
sidering the centroid results in the plane being forced to 
pass through the coordinate origin, i.e. a great circle fit is 
performed. This is done by most authors when calculating 
the moment of inertia tensor, but not here. We choose to 
seek planes without forcing them to go through the coor- 
dinate origin (the centre of the host galaxy), because this 
allows a consistency check: the constituent satellites must 
orbit within the host potential such that any disc made up 
of a virialised satellite population must pass near the origin. 
Test setups showed that the ALS algorithm is sensitive in 
deriving the correct distance of an artificial plane at an one- 
sigma level only. We calculate the centroid of the data-points 
and not the centre of mass of the satellites. Therefore, any 
found plane which has a distance to the host centre larger 
than about one disc-height may be interpreted as being un- 
physical. 

Since the ALS method is an unweighted fitting routine, 
distance uncertainties are accounted for by the applied error 
(AE) method: all satellites are randomly shifted along their 
line-of-sight with a normal distribution function as the prob- 
ability function of the magnitude of the shift. The variance of 
the normal distribution used for the shifting is derived from 
the distance uncertainty of the measurements. The random 
shifting is repeated a large number (10*) of times. For the 
analysis of the distribution of derived normals of the fit- 
ted planes the same analysis techniques as described below 
(§2.2) for the bootstrap re-sampling method are employed. 
Note that there are, however, schemes to introduce a weight- 
ing in ALS (e.g., Hartwick 2000; Zentner et al. 2005) but as 
a drawback this also influences the interpretation of the de- 
rived axis-ratios. 

A second method to fit planes is a weighted fit- 
ting routine based on the orthogonal distance regression 
(ODR) package provided by Netlib (http: //www .netlib . 
org/odrpack). An unweighted fit with the ODR method pro- 
vides the same result as the (much faster) ALS method. The 
weights are distance uncertainties in our application. This 
algorithm is a little bit subtle since it also estimates the 
errors of the fitted parameters. Even though the fitting pa- 
rameters converge, the error estimate may not, which can be 
understood as having weak constraints on the fitting param- 
eters. This only happened in special cases and never when 
applied to the satellites of the Milky Way or Andromeda 
directly. Further, the solution provided by ODR is strongly 
dependent on the weights used. The Cartesian variances cr^. 



ay, and derived from the distance uncertainty ar are not 
independent. Using the variances only leads to a different fit 
than using the full covariance matrix. Thus, we always use 
the full covariance matrix, accounting for correlations of the 
components. 

One measure of the planarity of the distribution is the 
flattening-parameter A/rcut as given in Kroupa et al. (2005). 
A is the root-mean-square height of the disc and rcut is 
the furthest distance to the Galactic Centre in the satel- 
lite sample. The alternative flattening-parameter A/rmcd as 
suggested by Zentner et al. (2005) is also a measure of the 
planarity, where Tmed is the median distance to the Galactic 
Centre of a sample. The formula to calculate A/rcut for an 
analytical r~'',q G R, distribution as given by Kang et al. 
(2005, their eqn. 1) is correct only for g < 3. The more 
general formula is given by: 



3-9 
3(5-q) 
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Z — q In |rcut I —In |ri 
3 



R, g/3,5 , 

g = 3, (1) 

g = 5 , 



where ri is the minimum radius of the distribution, 
liniri^o A always converges to for g ^ 3. Kroupa et al. 
(2005) derived a linear probability distribution p(r) oc , 
1.8 =^ p ^ 2.6 for the Milky Way (see also Koch & Grebel 
2006 for Andromeda), such that the spherical volume den- 
sity is psph (»-,■!?,<;/)) oc r-«, 3.8 ^ g sC 4.6. Kang et al. (2005) 
argued that the formally measured flattening A/rcut of a 
plane with inflnite number of particles following a power-law 
distribution decreases with power-law index g and converges 
to zero for g ^ 3. Indeed this always converges to zero for 
ri ^ 0, g ^ 3. A decreases with power-law index g, but more 
importantly, A is strongly dependent on the minimum ra- 
dius ri. This will also influence any test that is based on the 
measured height alone. Nevertheless, the argument by Kang 
et al. (2005) that a disc-like distribution may be mimicked if 
the satellite distribution is centrally concentrated with one 
or two outliers is valid (see also Zentner et al. 2005), re- 
quiring more robust statistical methods to be launched, this 
being one important aim of this study (§2.2 & §5). 

For the ALS method one can derive the axis-ratios c/a 
and b/a of the square-roots of the eigenvalues (ri ^ T2 ^ ra) 
of the moment of inertia tensor: c = y/ri, b = y/T2, a = y'rs. 
The values (a, b, c) are proportional to the root-mean-square 
(rms) deviation relative to the eigenvectors of To. In addi- 
tion we use the ratio c/b which indicates whether the tri- 
axial distribution is more oblate {c/b < b/a) or more pro- 
late {c/b > b/a). Note that this definition of 'triaxial more 
oblate' and 'triaxial more prolate' is different to the defi- 
nition based on the triaxiality parameter (e.g. Franx et al. 
1991). The ratio c/a is a better measurement than A/rcut 
and A / r^nod in terms of providing the ratio of the rms length 
in the direction of the smallest and largest extent, but can 
only be calculated for the ALS method. However, for a small 
sample, such as in the dwarf satellite application, quantities 
like A or the ratio c/a may not be a robust measure because 
of small number statistics. 
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2.2 The distribution of normal vectors of 
bootstrapped samples 

In previous works (Hartwick 2000; Kroupa et al. 2005; Kang 
et al. 2005; Libeskind et al. 2005; Zentner et al. 2005; Koch 
& Grebel 2006) the planarity of the satellite distribution was 
quantified based on the 'thickness' of the distribution in re- 
lation to some measure of the total spatial extent: A/vcut, 
A/r,ncd, or c/a. This is of course a basic requirement to call 
a distribution disc-like. However, if we think of a centrally 
concentrated distribution any small number of outliers will 
determine the orientation of the plane and we will always 
end up with a 'thin' disc (Kang et al. 2005; Zentner et al. 
2005). Consequently additional information about the ro- 
bustness of a disc-like distribution is needed to draw some 
statistically significant conclusions. This can be achieved us- 
ing a re-sampling technique. 

The bootstrapping method allows an estimate of the 
robustness of a disc-like distribution. The fitting to the re- 
sampled data is done using the unweighted ALS method. 
If the satellites are not distributed in a well-defined planar- 
like sheet, the normals of the fitted planes should show a 
large scatter. The amount of scatter of the directions of fit- 
ted normals to the bootstrapped samples is the quantity we 
use to determine the statistical significance of a plane-like 
distribution. To test the robustness of a best fitting plane to 
a set of data-points we quantify the spread of distributions 
of normals obtained from the large number of bootstrap- 
ping samples. A well-defined plane or disc of data points 
(satellites) will lead to a tight clustering of normals on the 
Galactic sky, while a weak disc will yield normals scattered 
over a large fraction of the sky. 

For detailed analysis of the bootstrapped data we fol- 
low methods described in Fisher, Lewis & Embleton (1987). 
The methods described below are only valid for a sample 
of unit vectors representing axial data, i.e. undirected data. 
For the present purpose the orientation of the normal of a 
plane is arbitrary, that is, the normal vector n represents 
the direction of an axis. A matrix M is defined as 



M 



^2 



yi 
m 



Z2 



(2) 



where {xi,yi,Zi) are the cartesian components of unit vec- 
tors fii, mis the number of unit vectors. An eigenvalue anal- 
ysis of the matrix 



(3) 



is performed. The eigenvector corresponding to the largest 
eigenvalue ts of T (n ^ r2 ^ ra) is the (estimated) principal 
axis ao of the input unit vectors, the normal vectors fii, and 
the principal axis corresponds approximately to the mean 
direction of the normal vectors. In our application, ao is the 
principal axis of the distribution of normals of fitted planes 
to the bootstrapped sample and we derive its direction on 
the Galactic sky. The shape parameter 7 is defined as 

_ In(r3/T2) 



ln(T2/ri) 

and the strength parameter ^ is defined as 
C = In (rs/n) 



(4) 



(5) 



(Fisher et al. 1987). These two quantities can be used to 
characterise a distribution of axial data on a sphere. 7 de- 
scribes the 'clusteriness' of the distribution. 7=1 indi- 
cates the transition between clustered (7 > 1) and girdled 
(7 < 1) distributions. is a continues parameter indicating 
the strength of concentration: the larger the more con- 
centrated a distribution is (being clustered or girdled), a 
uniform spherical distribution has C = 0. 

In addition the spherical standard distance 



sph 



Em [arccos(|ao • n,;j)] 



(6) 



is calculated where '•' denotes the scalar product of vectors. 
This is the analogue to the linear root-mean-square distance 
on the sphere. Since we deal with axial data we have to 
take the absolute value of the scalar product in Eq. (6). 
The spherical standard distance implies rotational symme- 
try but can be considered as an estimate of the upper limit 
of the opening angle of the sample of normal vectors for non 
rotational-symmetric distributions. 



3 DATA ANALYSIS 

For both large spirals of the Local Group, the Milky Way and 
Andromeda, satellite galaxies within the virial radius of the 
host galaxy are selected. This distance range is chosen since 
the satellite system within this radius is assumed to be viri- 
alised and the individual satellites are very likely bound to 
their hosts. Based on Local Group timing arguments (Kahn 
& Woltjer 1959) one may argue that the dark-matter haloes 
are more extended, that they even possibly overlap. This 
makes it difficult to associate some outer dwarf galaxies def- 
initely to one of the L* spirals or classify them as free floating 
objects. 



3.1 The Milky Way satellites 

For the Milky Way satellite system we use the data-set from 
Kroupa et al. (2005, their table 1), supplementing it with 
two newly discovered companions: the dwarf galaxies in Ursa 
major (UMa, Willman et al. 2005) at « 100 kpc and in 
Canes Venatici at a distance of 220l^g kpc (CVn, Zucker 
et al. 2006). Two other recently identified stellar systems, 
SDSS J1049-f5103 (Willman et al. 2005) and Bootes (Boo, 
Belokurov et al. 2006), are excluded from the present analy- 
sis because of their uncertain physical nature (globular clus- 
ter or dSph galaxy). We also exclude the Canis Major dwarf 
(CMa, Newberg et al. 2002; Martin et al. 2004) at 15 kpc 
from the Galactic Centre as its true nature is debated too 
(Moitinho et al. 2006). The latest discoveries of four further 
dSph satellites of the MW (Belokurov et al. 2006) were not 
included since the analysis for this paper was finished when 
the data were published. The positions of the Milky Way 
satellites are shown in Fig. 1 in an Aitoff projection. 



3.1.1 The KTB data-set 

First the data-set as used by Kroupa et al. (2005) is 
analysed. These are the most luminous satellites, compris- 
ing a complete census of satellites within the virial radius 
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Figure 1. An Aitoff projection of the positions for the innermost satelhtes within 254 kpe of the Milky Way as they would appear from 
the Galactic Centre (compare with Fig. 4 for M31). Also the position of Andromeda is shown. The positions of two recently discovered 
transitional objects are marked by open circles, and the positions of five further very recently discovered companions are marked by 
triangles. In this projection, the Sun is located at ^mw = 180°. The grey shaded region is the full sky-coverage region of the SDSS. 
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Figure 2. Positions of the MW satellites within the virial radius 
as seen from infinity. The top panels shows the results of the plane 
fitting using the unweighted ALS method and the bottom panels 
for the weighted ODR method when fitting through the complete 
sample of eleven brightest satellites marked by the filled dots (the 
KTB sample). In the left panels an edge-on view onto the fitted 
plane is shown, in the right panel a view rotated by 90° about the 
polar-axis of the Galaxy as indicated by the faint vertical dotted 
line. The MW is located at the origin and its disc is seen edge-on. 
The position of the UMa dwarf galaxy is marked with an open 
hexagon, the position of CVn with an open square. The shaded 
area shows the regions of fe = ±15° and b = ±30° which may be 
affected by obscuration through the MW disc. 



Table 1. The positions of the Milky Way satellite galaxies within 
the approximate virial radius of 254 kpc. In the first column we 
give a running number, in the second the name, in the third and 
forth longitude and latitude in galactocentric coordinates, and 
in the fifth the distance with 1-sigma errors from the Galactic 
Centre. In the sixth column the absolute luminosity in the V- 
band of the galaxies are given. 



No 


Name 


Imw 




''MW 


Lv 






n 


n 


[kpc] 




1 


Sgr 


8.8 


-21.5 


16.0 


± 2.0 




2 


LMC 


269.0 


-33.3 


50.2 


± 2.2 


2090(°) 


3 


SMC 


292.2 


-47.1 


56.9 


± 2.2 


575('=) 


4 


UMi 


114.5 


43.1 


68.1 


± 3.0 


0.29('>-) 


5 


Scu 


237.5 


-82.3 


79.2 


± 4.0 


2.15('') 


6 


Dra 


93.6 


34.7 


82.0 


± 6.0 


0.26('') 


7 


Sex 


237.1 


40.5 


89.2 


± 4.0 


0.50(=') 


8 


Car 


255.1 


-21.8 


102.7 


± 5.0 


0.43('') 


9 


For 


230.5 


-63.8 


140.1 


± 8.0 


15.5('') 


10 


LeoII 


216.5 


65.5 


207.7 


±12.0 


0.58^'') 


11 


Leol 


223.9 


48.1 


254.0 


±30.0 


4.79(a) 


12 


UMa 


162.0 


50.8 


104.9 


±20.0 


0.04('=) 


13 


CVn 


86.9 


80.2 


219.8 


±25.0 


0.12(d) 



References: Mateo (1998); Willman et al. (2005); ^"'^ Zucker 
et al. (2006); Lee & Kim (2000); van den Bergh (1999) 



of the Milky Way (rvir,Mw) s.nd brighter than Aftot.v ~ 
-8.8 mag, Lv = 2.6 x 10^ Lq (Mateo 1998). We fit a plane 
for the innermost eleven satellites (KTB data-set, without 
UMa and CVn) out to Leo I with a Galacto-centric distance 
of 254 kpc. For the ALS method, the direction of the nor- 
mal of the fitted plane is (Zmw = 157.3°, 6mw = —12.7°) 
and the distance of the plane from the Galactic Centre is 
Dp = 8.3 kpc, which appears to be quite large but is still 
well within the optical disc of the MW. An edge-on view and 
a view rotated by 90° about the polar axis of the Galaxy is 
shown in Fig. 2. The rms-height is A = 18.5 kpc, result- 
ing in flattening A/rcut ~ 0.07 and A/rmcci = 0.23. So the 
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distance of the fitted plane from the GC is a factor of two 
smaller than the rms-height of the plane. Note that Kroupa 
et al. (2005) found A/rcut = 0.10, and Zentner et al. (2005) 
found A/rniod ~ 0.3 for the same data. The larger values 
than derived here are caused by the suboptimal fitting rou- 
tines used there. The derived axis ratios are c/a = 0.18 and 
b/a = 0.53, resulting in c/b = 0.34 < b/a, a tria^xial, highly 
oblate distribution of satellites. Libeskind et al. (2005) gave 
values of c/a ~ 0.3, b/a ~ 0.5 for the MW in their figure 3, 
which would mean c/b — 0.6 > b/a, i.e. a triaxial, slightly 
prolate distribution. 

Next, the applied error (AE) method is applied, ran- 
domly shifting the position of all satellites along their line-of- 
sight vector with a normal distribution function as the prob- 
ability function of the magnitude of the shift, repeating this 
10^ times. The principal axis of the resulting distribution 
of the normals is located at {Imw = 157.4°, 6mw ~ —12.6°) 
with a spherical standard distance Asph = 1.2°. The derived 
principal axis is in good agreement with the single fit above 
and the scatter as quantified by A^ph is remarkably small 
about the principal axis. 

Using the ODR method the pole is located at {Imw = 
158.2°, 6mw = -29.0°) with a distance Dp = 3.4 kpc of 
the plane from the GC. The rms-height is A = 32.6 kpc, 
resulting in flattening A/rcut = 0.13 and A/rmcd = 0.41 
The longitude of the derive pole is very similar for the ALS 
and the ODR method, while the latitude deviate by ~ 16°. 
This can be understood by looking edge-on onto the fit- 
ted plane (Fig. 2, left panels). Since we are basically sit- 
ting in the plane of the satellites, most of the distance un- 
certainties, which are considered in the ODR method, are 
along the radius of the plane. As can be seen in Fig. 2, for 
the distant satellites the components of the distance uncer- 
tainties along the polar-axis (the ordinate) are the largest. 
This forces the ODR algorithm to weight positions along the 
polar-axis down which results in a different latitude of the 
pole while the longitude of the fitted normal is not affected. 
This outcome is an indication that, while the ODR method 
is robust against single outliers, it can be biased strongly by 
a systematic alignment of the provided distance uncertain- 
ties. Later we show that this aspect affects the fitting for 
the Andromeda satellites even more. 

To investigate the robustness of the disc-like feature, 
10 000 bootstrap re-samplings for the eleven innermost satel- 
lites of the Milky Way (^A^tot = 352 155 for n = 11, see 
Eq. A2) are performed with the ALS method. In Fig. 3 a 
smoothed {Imw, 6mw) scatter plot of the locations of the 
normals of the bootstrap samples is shown in grey-scale, 
white corresponding to zero density. The smoothing ker- 
nel is a Fisher function with smoothing parameter k = 100 
(Fisher et al. 1987)^. In addition the contour lines for the 
density estimate of 0.5, 0.75 and 0.9 are plotted. The plot is 
centred on the principal axis of the distribution which is lo- 
cated at (Zmw ~ 158.2°, 6mw = —11.9°) and thus is in good 
agreement with the results of the single fit above. The shape 
parameter 7 = 3.8 and the strength parameter ^ = 4.3 show 
that the distribution is strongly clustered around its princi- 



The Fisher function is the equivalent of the Gaussian on a 
sphere. The larger the smoothing parameter k is chosen the nar- 
rower the smoothing kernel gets. 




Figure 3. A smoothed (^mw, ^mw) scatter plot of the distribu- 
tion of normals of planes fitted to each of 10 000 bootstrap sam- 
ples for the eleven innermost MW satellites (without the recent 
additions). A spherical density estimate is shown using a Fisher 
density for the kernel function (k = 100, Fisher et al. 1987) plot- 
ted in grey-scales, white corresponds to 0. In addition the contour 
linos for the density estimate of 0.5, 0.75 and 0.9 in units of the 
central density are plotted. The principal axis of the distribution 
is located at {Imw = 158.2°, 6mw = -11.9°). 



pal axis, and well outside the distribution of 7, 1^ values for 
an intrinsically isotropic distribution of satellites (Fig. 11). 
The resulting spherical standard distance is Agph = 13.0°. 

The spherical standard distance Asph found for the AE 
test is only a tenth of that found with bootstrapping. This 
clearly shows that the systematic error caused by the dis- 
tance uncertainties of the Milky Way satellites is signifi- 
cantly smaller than the intrinsic scatter determined with the 
bootstrapping. So the disc-like feature is not affected much 
by the distance uncertainties of the MW satellites. 



3.1.2 Including the newly discovered satellites 

We repeat the same analysis, now consecutively including 
the UMa and CVn dwarf galaxies. The results from the sin- 
gle flts, the AE analysis and the bootstrap analysis are given 
in Tables 2-4. For the flts the locations of the poles are only 
marginally affected since both dSphs are located close to the 
former fitted disc (see Fig. 2). With a distance of ~ 220 kpc 
CVn is in fact the second furthest satellite galaxy in our 
sample. For the unweighted ALS method the orientation of 
the fitted disc is mostly determined by the outer satellites, 
nevertheless the orientation is not affected much. 

Including both new satellite galaxies, the spherical stan- 
dard distance Asph found for the AE test remains an order 
of magnitude smaller than for the bootstrapping, showing 
that the distance uncertainties do not systematically affect 
our results. 
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Table 2. Results from the single plane fits for the innermost satellite galaxies of the Milky 
Way within the approximate virial radius for the two fitting methods: ALS and ODR (see 
Text). Results for the KTB data-sot, and for the data-sets including the UMa and the CVn 
dwarf galaxies are given. For Andromeda the results for both data-sets, MI and KG, as well 
as for a morphologically (mss8, §3.3) and a kincmatically motivated (kss8, §4.3) subsamplc 
of eight satellites are tabulated. Results are given in Galacto-centric and Andromeda-centric 
coordinates, respectively. 



method 


'mw [°] 


''MW [°] 


Dp [kpc] 


A [kpc] 


A/rcut 


A/rmcd 


c/a 


b/a 








Milky Way: KTB 










ALS 


157.3 


-12.7 


8.3 


18.5 


0.07 


0.23 


0.18 


0.58 


ODR 


158.2 


-29.0 


3.4 


32.6 


0.13 


0.41 










Milky Way: KTB -f UMa 








ALS 


160.5 


-14.6 


12.5 


20.3 


0.08 


0.25 


0.21 


0.54 


ODR 


158.1 


-29.1 


3.4 


31.7 


0.12 


0.39 










Milky Way: KTB + UMa -f CVn 








ALS 


153.8 


-10.2 


7.8 


22.8 


0.09 


0.27 


0.22 


0.55 


ODR 


157.4 


-29.2 


2.9 


40.9 


0.16 


0.48 








'm31 [°] 


^M31 [°] 


Dp [kpc] 


A [kpc] 


A/rcut 


A/rjnod 


c/a 


h/a 








Andromeda, 


Ml-data 










ALS 


73.4 


-31.5 


1.0 


45.9 


0.17 


0.42 


0.36 


0.46 


ODR 


23.8 


-12.5 


45.0 


54.4 


0.20 


0.50 






ALS mss8 


177.0 


-24.1 


34.9 


29.2 


0.11 


0.39 


0.27 


0.67 


ODR mssS 


177.9 


-20.4 


37.7 


29.5 


0.11 


0.39 






ALS kss8 


69.9 


-35.2 


1.8 


16.5 


0.06 


0.11 


0.12 


0.50 


ODR kssS 


57.9 


-28.6 


12.8 
Andromeda, 


22.7 
KG-data 


0.08 


0.16 






ALS 


83.5 


-31.0 


7.5 


46.1 


0.16 


0.46 


0.41 


0.68 


ODR 


27.6 


-31.1 


32.5 


68.2 


0.24 


0.68 






ALS mss8 


168.0 


-26.7 


1.6 


9.4 


0.04 


0.14 


0.09 


0.68 


ODR mssB 


168.5 


-29.4 


1.9 


10.2 


0.04 


0.15 






ALS kss8 


73.6 


-35.0 


3.9 


17.9 


0.06 


0.18 


0.15 


0.71 


ODR kssS 


52.9 


-35.5 


13.7 


31.3 


0.11 


0.31 







Table 3. Results from the applied error test for satellite galaxies 
of the Milky Way and Andromeda. The longitudes and latitudes 
of the derived principal axis and the spherical standard distance 
of the distributions for the data-set used by Kroupa et al. are 
listed. The recently discovered dSph galaxies in Ursa Major and 
Canes Venatici are included. The results for the two data-sets 
used for Andromeda, as well as for two subsamples (§3.3 & §4.3) 
of M31 satellites are given . 



Data-set 


'mw [°] 


hmw [°] 


Asph [°] 


Milky Way KTB 


157.4 


-12.6 


1.2 


-l-UMa 


160.6 


-14.6 


1.5 


-l-CVn 


153.9 


-10.2 


1.5 




Imzi [°] 


6m31 [°] 




Andromeda Ml 


75.1 


-31.7 


13.2 


mss8 


176.6 


-24.9 


12.5 


kssS 


70.5 


-35.2 


2.4 


Andromeda KG 


82.9 


-31.1 


8.0 


mss8 


165.2 


-30.9 


21.5 


kssS 


74.2 


-35.2 


1.6 



3.2 The Andromeda satellites 

In order to study the three-dimensional distribution of the 
satellite system of Andromeda it is most convenient to 
transform their position vectors relative to the observer 
into an Andromeda-centric coordinate system (see also Mc- 
Connachie & Irwin 2006b). A detailed description of the 
transformation is given in Appendix B in a general way. Two 
different data-sets for the distances of Andromeda and its 
satellites were incorporated: the first data-set as published 
by McConnachie & Irwin (2006b, MI data-set, see their table 



Table 4. Results from the bootstrap re-sampling method for the 
satellite galaxies of the Milky Way and Andromeda. The same 
quantities as in Table 3 are provided and in addition in columns 
five and six the shape parameter 7 and the strength parameter ^ 
of the distribution of normals of the fitted planes are tabulated. 



Data-set 


'mw [°] 


^IMW [°] 


Asph [°] 


7 


c 


MW KTB 


158.2 


-11.9 


13.0 


3.8 


4.3 


-l-UMa 


161.8 


-14.9 


12.0 


4.0 


4.4 


-l-CVn 


156.7 


-10.6 


12.4 


2.8 


4.6 




'm31 [°] 


bm-ii [°] 








M31 MI 


75.5 


-31.9 


38.6 


0.6 


3.1 


mss8 


178.3 


-28.5 


32.7 


1.1 


2.7 


kssS 


69.5 


-34.2 


9.8 


5.9 


4.7 


M31 KG 


83.1 


-30.0 


27.9 


2.9 


2.6 


mss8 


167.1 


-29.6 


11.8 


5.3 


1.9 


kssS 


72.4 


-33.8 


11.5 


14.9 


4.2 



1) where most of the distances were derived using the tip of 
the red giant branch method using ground based telescopes 
(McConnachie et al. 2005). The other data-set as given by 
Koch & Grebel (2006, KG data-set, see their table 1): they 
compiled a list of HST-based distance measurements. The 
data is given in Table 5 for both data-sets in Andromeda- 
centric coordinates. 

Fig. 4 shows an Aitoff projection of the satellite distri- 
bution on the Andromeda sky. The error bars due to the 
combined uncertainties in the distance measurement of the 
satellites and Andromeda are shown for both data-sets: in 
grey with thick lines for the KG data-set and in black for 
the MI data-set. Note the voids in the regions Im-h < 180°, 
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-90° 



Figure 4. The AitofT projection of the M31 sateUites in the Andromeda-centric coordinate system for two data-sets: in grey with thick 
Hnes for the data given in Koch & Grebel (2006) and in black for the data given in McConnachic & Irwin (2006b). The error bars due 
to the combined uncertainties in the distance measurement of the sateUites and M31 are shown. The position of the MW is also marked. 



feM3i < 0° and /m3i > 180°, feM3i > 0°. As shown by Mc- 
Connachie & Irwin (2006b) these regions are only marginally 
obscured by the MW disc (see their figure 1, the region of 
maximum obscuration is near IC 10). 

Similar to the case of the MW, more faint dwarf galax- 
ies probably remain to be found for M31 within the next 
few years. One of these discoveries was recently reported by 
Zucker et al. (2006) : the dSph And X is comparable in lumi- 
nosity to And IX, however the distance determination was 
difficult. Zucker et al. (2006) gave a distance of 667±30kpc 
to 738±35 kpc. Another three satellites. And XI - And XIII, 
were reported in a very recent paper (Martin et al. 2006). 
No distances could be determined for the individual satel- 
lites, but combining their colour-magnitude diagrams and 
assuming all to have the same distance, Martin et al. (2006) 
derived a combined distance of 740 - 955 kpc. Given these 
uncertainties we do not include the four new satellites in our 
analysis but discuss them later. 

The innermost twelve satellites of Andromeda without 
M33 lie within « 269 kpc which is the approximate virial 
radius of Andromeda (rvir,M3i). We assume LGS 3 to be the 
twelfth satellite, although for the KG data-set And VI is ac- 
tually closer to the centre. This will be addressed later (this 
section and §4.2). The results are given in Table 2. For the 
Ml-data the distribution is triaxial and more prolate, while 
for the KG data-set the distribution is found to be triaxial 
and more oblate, which is reflected in the further analysis, 
too. Fig. 5 shows an edge-on view of the fitted planes and a 
view rotated by 90° about the polar axis of Andromeda for 
both data-sets. As can be seen the recently discovered dSph 
And X that was not incorporated in the fitting (marked by 
the open circle near the centre of the plots) is located close to 
the fitted plane. If we just recalculate the rms-heights of the 
fitted disc, now including And X, it increases only slightly 
to A = 46.6 kpc for the MI data-set, A = 47.2 kpc for the 
KG data-set. Interestingly, also the three very recently dis- 
covered satellites And XI - And XIII (Martin et al. 2006) 
all lie very close (« 10 kpc, [« 30 kpc]) to the fitted disc 
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Figure 5. The Andromeda satellite system as seen from infinity. 
An edge-on view (left panel) and a view rotated by 90° about the 
polar axis (right panel) of the fitted planes derived using the ALS 
method arc shown as in Fig. 2. The top panels are the results 
for the KG-data and the bottom for the Ml-data. The different 
symbols mark different morphological types: dots mark dSphs, 
squares cE/dEs, and triangles dIrr and dIrr/dSph galaxies. Those 
satellite galaxies marked with filled symbols were incorporated in 
the fitting routine while those with open symbols were left out. 
M33, which was also not incorporated in the fitting, is marked by 
the pentagram. Satellites encircled by a light grey circle were later 
excluded from fitting to a kinematically motivated subsample, 
see §4.3. The error bars given are derived from the line-of-sight 
distance uncertainties of the satellites only. 
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Table 5. Basic parameters of the Andromeda satellites: the first four columns contain the a running number, the name, the 
morphological type and the absolute luminosity in the V-band of the satellite, columns five to seven list the positions in 
Andromeda-centric coordinates (§B1) for the first data-set (McConnachie & Irwin 2006b), and the eighth and ninth columns 
list the radial and perpendicular components of the measured line-of-sight velocity relative to Andromeda (§4.1). In the columns 
ten to fourteen the same data are provided for the second data-set (Koch & Grcbcl 2006). The asterisks mark satellites for 
which the possible poles of the angular momentum vector can be restricted to an arc of 180° (§4.1). 

Ml-data KG-data 



No 


Name 


Type 


Lv 


'm31 


^M31 


''M31 


Vt 




'm31 


^M31 




ft 










[lO^Lo] 


[°] 


n 


[kpc] 


[kms -*-] 


[°] 


[°] 


[kpc] 


[kms-i] 


1 


M32 


cE 


383'-='' 


278.5 


-35.7 


6 


96.0 


-1.8* 


242.3 


-37.9 


6 


83.2 


-47.9 


2 


NGC 205 


dE 


366(=') 


0.9 


25.1 


40 


15.0 


55.1 


0.5 


21.3 


58 


11.3 


55.9 


3 


And IX 


dSph 


0.17 f'^) 


118.3 


-11.1 


42 


68.6 


-51.6 


65.4 


-0.2 


40 


82.2 


24.7 


4 


And I 


dSph 


4,37{b) 


220.1 


-42.1 


59 


71.4 


38.9 


306.2 


-37.8 


48 


66.9 


-46.3 


5 


And III 


dSph 


0.58*'') 


241.1 


-24.1 


76 


61.0 


0.6 


260.0 


-22.8 


68 


58.4 


-17.6 


6 


And V 


dSph 


0.58('') 


99.6 


0.2 


110 


106.2 


-38.9* 


76.0 


5.3 


117 


81.0 


-78.9 


7 


NGC 147 


dE 


13lW 


155.8 


16.1 


145 


46.1 


-102.5 


122.0 


42.3 


101 


100.3 


-50.7* 


8 


And II 


dSph 


9.12('') 


188.7 


-50.9 


185 


47.1 


-111.3 


193.9 


-60.9 


160 


66.4 


-101.0 


9 


NGC 185 


dE 


125(=') 


162.4 


3.9 


190 


17.5 


-102.1 


160.8 


5.5 


175 


21.5 


-101.4 


10 


And VII 


dE 


27.5('') 


170.1 


63.3 


219 


14.4 


-84.5* 


169.2 


65.6 


216 


11.0 


-85.0 


11 


IC 10 


dirr 


160(=') 


103.2 


45.2 


260 


26.3 


-108.2 


140.6 


27.2 


255 


83.3 


-73.9 


12 


LGS 3 


dl/dS 


1.33 


264.0 


-53.2 


269 


8.8 


-102.1 


220.4 


-45.9 


284 


41.3 


-93.8 


13 


And VI 


dSph 


3.3l('') 


259.8 


-3.5 


269 


53.0 


-112.1* 


260.6 


-3.3 


266 


51.4 


-112.8* 


14 


Peg DIG 


dl/dS 


12.0(^) 


270.1 


-1.0 


474 


145.5 


-93.0 


252.9 


-4.7 


410 


110.7 


-132.6* 


15 


IC 1613 


dIrr 


63.6*^) 


243.1 


-46.7 


511 


39.2 


-187.3 


244.8 


-46.9 


505 


43.2 


-186.4 


16 


M33 


Sc 


3020(<^) 


340.2 


-77.6 


207 


126.3 


-47.6* 


350.0 


-65.2 


221 


133.7 


-18.5 



References: Mateo (1998); C') McConnachie & Irwin (2006a); Zucker et al. (2004); (d) van den Bcrgh (1999) 



when using a mean combined distance of 847.5 kpc albeit 
with large uncertainties. 

For the applied error method (see Table 3) the distance 
of Andromeda is not varied. This would only affect the dis- 
tance of the plane from the origin. Compared to the spher- 
ical standard distance derived for the MW, it is a factor of 
ten larger for the M31 system, which is not surprising given 
the much larger uncertainties in the relative distances M31- 
satellite. The principal axes come out in good agreement 
with the poles of the single fits. 

However, the ODR method yields poles that are far 
away from the poles found with the ALS method (Table 2). 
The difference is totally dominated by the large distance 
uncertainties, which are here not aligned within the fitted 
plane as in the case of the Milky Way. Instead they are sys- 
tematically aligned along the LOS from the MW to M31. 
The components in the direction to the MW are weighted 
down and the fitted plane appears to be nearly perpendic- 
ular to the direction from the MW to M31, a result of the 
systematic dependencies of the covariances. 

McConnachie et al. (2005) showed that the satellite dis- 
tribution is significantly offset towards the direction of the 
Milky Way (also visible in Fig. 5, right panels). The offset is 
reflected by the large distance from the centre of M31 along 
the direction of the normal for the fitted planes when using 
the ODR method due to the nearly face-on orientation of the 
fitted disc. In contrast, the disc-like distribution found with 
the ALS method is more edge-on, i.e. the systematic offset 
as identified by McConnachie et al. is within the plane. 

As for the MW 10 000 bootstrap re-samplings for both 
data-sets of M31 are performed to test the robustness of 
the plane. The principal axis of the distribution is (/m3i = 
75.5°, feM3i = -31.9°) [(83.1°, -30.0°)], being in good 
agreement with the original fit. We derive a shape parame- 
ter 7 = 0.6 [2.9] and a strength parameter = 3.1 [2.6], the 



spherical standard distance is A^ph — 38.6° [27.9°] (num- 
bers for the MI[KG] data-sets). While for the KG data-set 
the distribution of the directions of fitted normals for the 
bootstrapped sample is found to be clustered (7 > 1), for 
the Ml-data it is found to be a girdled distribution (7 < 1). 
Fig. 6 shows a smoothed (Zm3i, feM3i) scatter plot of the dis- 
tribution of the fitted normals for the MI data-set: there is 
a distinct peak about the principal axis and a second, very 
weak over-density can be seen nearly 90° off, being the ori- 
gin of the girdled distribution. The KG data-set does not 
show a secondary maximum. 

And VI has approximately the same distance from M31 
as LGS 3, both lying close to the approximate virial radius. 
Including And VI in the fitting routine for the KG-data, 
where it is actually closer to the centre of M31 than LGS 3, 
dramatically changes the picture. The pole of the fitted nor- 
mals is clearly offset from the fits without And VI, the dis- 
tance of the fitted plane is significantly offset from the centre 
of Andromeda, and the axis ratios do change significantly: 
b/a = 0.63 and c/a = 0.57. More importantly the clustering 
found for the bootstrapping without And VI for the KG-data 
disappears completely. If the satellite was within a planar- 
like distribution, the bootstrapped distribution should be- 
come similarly or even more tightly concentrated as it did 
when including the UMa dwarf galaxy for the Milky Way. 
Instead, it gets very weak: 7 = 0.3 and = 1.7. Therefore 
we treat And VI as an outlier. 

In contrast to the Milky Way satellite system, for An- 
dromeda the spherical standard distance derived with the 
applied error method is of the same order as for the boot- 
strap method which is a result of the large distance uncer- 
tainties for M31 and its satellites. So the results may well 
be affected by the still too large distance uncertainties for 
M31 and its satellite galaxies. 
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Figure 6. A smoothed (Zm31, ^msi) scatter plot of 10 000 boot- 
strap samples for the innermost twelve Andromeda satellites for 
the MI data-set ranging out to LGS 3 as Fig. 3. The principal axis 
of the distribution is located at (Zm31 = 75.5°, 6m31 = —31.9°). 
The plot is shown 30° off-centre from the principal axis. An ad- 
ditional contour line for the density estimate of 0.25 indicates a 
very weak secondary maximum. 
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Figure 7. An edge-on view (left panels) and a view rotated by 90° 
about the polar axis (right panels) of the fitted plane (using ALS) 
for the morphologically motivated sub-subsample as proposed by 
Koch & Grebel. In the top panels we show the results for the MI 
data-set and in the bottom for KG data-set. Symbols are chosen 
as in Fig. 5. Satellites not incorporated in the fitting are marked 
with open symbols. The error bars given are derived from the 
distance uncertainties of the satellites only. 



3.3 Morphological subsample of Andromeda 
satellites 

Koch & Grebel (2006) found a very pronounced polar disc- 
like feature for a morphologically motivated subsample of 
early-type dwarf galaxies. Their procedure was as follows: 
they first fitted a plane to all seven dSph galaxies in their 
data-set and then excluded And II (at a distance of tmsi = 
160 kpc < rvir,M3i) as an outlier because of its large dis- 
tance to the fitted plane. This disc-like feature of six dSphs 
was found to be highly statistically significant. Next they in- 
cluded all other morphologically similar galaxies, three dEs 
and one cE (M32), again finding a disc-like feature with high 
statistical significance. In a last step they excluded two of 
the three dEs, but included one transitional type object, the 
dIrr/dSph Peg DIG at a distance of 410 kpc > rvir.Msi, be- 
cause of its close proximity to the disc-like feature found 
before'^ . 

Indeed we confirm an amazingly thin configuration 
when using this sub-subsample consisting of M32, And I, 
And III, NGC 147, And V, And VII, And VI, And IX, and 
Peg DIG. In Tables 2 - 4 we refer to the sub-subsample 
without Peg DIG as mss8 (morphological subsample of eight 
satellites). We exclude Peg DIG because it is well outside the 
approximate virial radius of M31. The fitted configuration 
is shown in Fig. 7. For the KG data-set the pole of the fit- 
ted plane is (Zm3i = 168.0°, 6m3i = —26.7°), with a distance 



^ However, also note the different scaling of the axes in their 
figure 3 which makes the distributions appear more planar-like 
than they truly are. 



from the centre Dp = 1.6 kpc, a rms-height A = 9.4 kpc, and 
with axis ratios c/a = 0.09 and b/a = 0.68. For the weighted 
ODR method the results are very similar (see Table 2) . M33 
is located very close to this fitted plane. Including Peg DIG 
for completeness results in (Zmsi = 163.0°, feM3i ~ —27.3°), 
with a distance from the centre Dp — 1.2 kpc, a rms- 
height A = 13.1 kpc, and with axis ratios c/a = 0.09 and 
b/a = 0.45. But this distribution (6m3i = —27.3°) is not as 
polar aligned as claimed by Koch & Grebel (2006) due to 
the incorrect transformation to the Andromeda-centric co- 
ordinate system used by them (§ Bl). Also note that from 
the Milky Way we are basically looking face-on onto this 
fitted plane, which is an important clue as we show later. 

Using the MI data-set without Peg DIG (MI mss8) leads 
to (/m3i = 177.0°, 6m3i = -24.1°), Dp = 34.9 kpc, and 
A = 29.2 kpc; including Peg DIG, (Zm3i = 182.4°, 6m3i = 
-23.2°), Dp = 35.1 kpc, and A = 29.1 kpc. For this data- 
set the fitted plane is not as thin as for the KG data, and 
the offset from the centre of M31 is remarkably larger than 
the rms-height of the fitted plane. As can be seen in Fig. 7, 
lower panel, there is now another dE (NGC 147, the filled 
square to the right of the plane) remarkably offset from the 
fitted plane. 

Further insight comes from the AE test (Table 3). When 
appUed to the morphological sub-subsample for the derived 
spherical standard distance is Asph = 21.5°. This is a very 
large uncertainty in the location of the poles of the fitted 
normals, a factor three larger than for the full set of twelve 
satellites within the approximate virial radius used before. 
For the MI data-set, Agph = 12.5° is of the same order as 
for the full data-set. 
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To test the robustness of the results a bootstrap analysis 
for the sub-subsample without Peg DIG is performed, now 
using 5 000 re-samplings accounting for the smaller number 
of possible distinct bootstrap samples ('^A'tot = 6 231 for 
n = 8). The results are given in Table 4. The principal axes 
are in agreement with the single fits. The distribution of di- 
rection of bootstrapped normals is marginally clustered and 
concentrated. However the spherical standard distance for 
the bootstrapped sample is remarkably smaller than for the 
AE test. This indicates that the systematic uncertainties of 
the distances are larger than the intrinsic scatter of the fitted 
disc for the KG data. For the MI data-set the bootstrapped 
distribution is not found to be clustered but of transitional 
type. The spherical standard distance is significantly larger 
than for the KG data-set. 

The recently discovered dSph And X is also found to 
be off the disc-like structure of the mss8 subsample for the 
KG data-set. For a heliocentric distance of 702.5 kpc its 
distance from the fitted disc is « 87 kpc (> 9(t, And II, 
excluded as an outlier by Koch & Grebel, is « 127 kpc 
away). To be within ±10 kpc from the disc, And X would 
have to be at a heliocentric distance of « 786 - 808 kpc. 
For the MI data-set And X's distance to the fitted disc is 
« 32.7 kpc. And X is located on the near side of M31 to 
the MW, thus adding to the systematic offset of the M31 
satellite system towards the barycentre of the Local Group 
(McConnachie & Irwin 2006b). 

We thus find that the apparent disc-like configuration of 
the dSph/dE satellites sub-subsample for M31 is present for 
the KG data-set only and can not be reproduced using the 
MI data-set. The nearly face-on alignment relative to the 
Milky Way results in distance uncertainties basically per- 
pendicular to the fitted plane (Fig. 7). The thin configura- 
tion disappears when shifting the satellites along their line- 
of-sight in accord with the distance uncertainties as done 
in the AE test. Comparing the results for the AE test and 
the bootstrapping suggests that the systematic uncertainties 
caused by the distance-measurement errors are larger than 
the intrinsic scatter of the distribution. Thus the thin disc- 
like configuration found may be just a chance alignment for 
the KG data-set, but its existence is not completely ruled 
out. 



4 KINEMATICS OF THE M31 SATELLITES 
4.1 Restricted polar paths 

Based on geometrical arguments, Lynden-Bell & Lynden- 
Bell (1995) constructed a diagram of so-called 'polar paths' 
for the Milky Way satellite galaxies and globular clusters. 
The direction of every possible pole of a satellite orbit must 
be located at a right angle to its direction from the Galac- 
tic Centre. Regions in the sky where three or more polar 
paths (nearly) intersect define the poles of possible streams 
of satellites. Palma et al. (2002) used proper motion mea- 
surements to restrict the possible orbits, which resulted in 
arc segments of possible poles rather than great circles. 
McConnachie & Irwin (2006b) used the available three- 
dimensional spatial data for the M31 satellites to construct 
polar paths for Andromeda and identified multiple possi- 
ble stream candidates of satellite galaxies. We note that a 




Figure 8. Geometrical illustration how to calculate radial and 
perpendicular velocity components of the measured line-of-sight 
(LOS) velocity of an Andromeda satellite. In the top panel we 
show the measured vectors in the rest-frame of the Sun. The ve- 
locity perpendicular to the line-of-sight, i't,L0Si is unknown, i.e. 
there is a family of possible total velocity vectors, indicated by 
dotted arrows. In the bottom panel we show the velocities cor- 
rected for the LOS velocity of Andromeda. The corrected LOS 
velocity of a satellite can be split into a radial component Vr and 
a perpendicular component vt relative to M31. 



a disc-of-satellites is rotationally supported if the angular 
momenta of its constituent satellites are aligned. Since in 
a non-spherical halo the angular momentum vectors would 
precess, the mutual distances of the directions of the angular 
momentum vectors will grow in the course of time. 

The observer's view of the Andromeda satellite system 
is very different from that of the Milky Way. While the MW 
satellites are basically seen from the Galactic Centre and 
only the radial velocity component is available, the M31 
satellite system is almost seen from infinity. Together with 
the spatial information relative to M31 it is therefore possi- 
ble to decompose the observed line-of-sight (LOS) velocities 
into radial and perpendicular components relative to M31. 
These are, however, only lower limits of their total velocities 
since only the LOS velocities can be measured. 

Fig. 8 illustrates how the measured line-of-sight velocity 
of an M31 satellite breaks up into radial and perpendicular 
components relative to M31. The top panel shows the mea- 
sured heliocentric velocity vectors ulos and ?;los(M31). In 
the bottom panel the same velocity vectors as in the top 
panel are shown but now corrected for the LOS velocity of 
M31, i.e., in the rest-frame of M31 assuming no consider- 
able perpendicular velocity (proper motion) for Andromeda 
(Kahn & Woltjer 1959; Einasto & Lynden-Bell 1982), aL 
though there are arguments which allow for a significant 
proper motion component of M31 (Loeb et al. 2005). The 
velocity component dt,los perpendicular to the line-of-sight 
of a satellite galaxy is unknown. Therefore a whole family 



12 M. Metz, P. Kroupa and H. Jerjen 



of velocity vectors is possible for the satellite, illustrated by 
dotted vectors. 

The transformation of velocity vectors into the 
Andromeda-centric rest-frame is calculated by taking the 
time derivative of Eq. (B4) : 



VM31 = Rm3i (v - t;(M3l)) 



(7) 



Upon setting i;(M3l) = i'los(M31) the line-of-sight velocity 
vectors of the satellites are calculated in the Andromeda- 
centric coordinate system. The LOS velocity vector of the 
satellite can now be split into a radial and a perpendicular 
component, Vr and vt, respectively. The values of the ra- 
dial and perpendicular velocity components relative to An- 
dromeda are given in Table 5 for the M31 satellites, negative 
radial velocity meaning that the component is pointing to- 
wards M31, positive that it is pointing away. 

This information is used next to set some limits on the 
possible poles of the orbits of the M31 satellites. In the ex- 
ample shown in Fig. 8 without an additional velocity com- 
ponent perpendicular to the line-of-sight, the sense of rota- 
tion about M31 would be counterclockwise (within the MW- 
M31-satellite plane). Only a large enough velocity compo- 
nent vi perpendicular to the LOS and within the MW-M31- 
satellite plane can reverse the sense of rotation. Any compo- 
nent ut,los perpendicular to the MW-M31-satellite plane 
cannot reverse the sense of rotation, but can only displace 
the direction of the angular momentum vector along the po- 
lar path by a maximum ±90° from the direction derived 
assuming no perpendicular velocity. 

We have to make some additional assumptions about 
the maximum allowed velocities of the satellites. We assume 
the satellites are bound to an isothermal dark-matter halo 
of M31 with circular velocity Vc ~ 250km s^^, truncated at 
Vc = 250 kpc. The escape velocity is given by 



Ve{r) = 



\/2 Uc VI + In fc — In r 



r < rc 
r > rc 



(8) 



A velocity component perpendicular to the LOS in the 
plane of the MW, M31, and the satellite can either con- 
tribute to the perpendicular velocity (giving a total velocity 
vector V2) or counteract (i^i). We calculate the maximum 
possible velocity vectors v-i and V2 in the M31 rest-frame 
with absolute values of the velocities = \v2\ ~ Ve{rs) 
(Eqn. 8), where Vs is the distance of a satellite from the cen- 
tre of M31. The direction of the orbit of a satellite can now 
be restricted if the perpendicular components of i;i and V2 
both point in the same direction, i.e., if 

Vi ■ Vt V2- Vt 



\Vl\\Vt\ \V2 Vt 



= +1 



(9) 



where '•' denotes the scalar product of vectors. This allows 
restriction of the polar paths to at least arcs of 180° which 
we call "restricted polar paths" (RPPs), i.e., the orientation 
of the orbit about M31 can be inferred. 



4.2 Application of restricted polar paths 

RPPs are calculated using heliocentric radial velocities as 
given in McConnachie & Irwin (2006b, their table 1) for 
both data-sets. In Table 5 we mark the data for those satel- 
lites with an asterisk for which we can restrict the poles by 



the RPP method described above. Only for And VI the poles 
can be restricted for both data-sets. The polar paths and the 
resulting RPPs are plotted in Fig. 9 for the MI data and in 
Fig. 10 for the KG data. Three projections for each data-set 
are displayed: an Aitoff-projection shows the overall global 
distribution of the polar paths, while it is difficult to dis- 
tinguish the different tracks near the poles. For this reason 
two Lambert-projections are plotted in addition, in the top- 
left for the northern hemisphere, and in the bottom-right 
for the southern hemisphere. Data for the different satel- 
lites are plotted in different colours (see online material for 
coloured version). The polar paths for those satellites which 
were used for fitting the plane in Sect. 3.2 are plotted with 
solid lines (satellites 1 - 12 in Table 5), for satellites not 
used in the fitting with dashed lines. The position of the 
angular momentum vector of M33 as derived from the mea- 
sured radial velocity and proper motion (Brunthaler et al. 
2005) is marked by a triangle assuming no proper motion 
for Andromeda. A considerable proper motion component 
of M31 would shift the direction of the angular momentum 
vector of M33 along its polar path in Figs. 9 and 10 since 
the measured velocity of M33 is transformed to the rest- 
frame of M31. In addition we plot loops with distances of 
15° and 30° from the pole of the fitted plane using the ALS 
method for all satellites within the approximate virial radius 
of M31 (§3.2). Since the direction of the normal is arbitrary 
for the fitted plane, two regions appear on opposite sides of 
the sphere. 

There is a strong clustering of intersection points of 
polar paths near the pole of the fitted plane, both in the 
northern and southern hemisphere. Those satellites whose 
polar paths intersect there are candidates for kinematic 
streams. In particular these are: NGC 205, And I, NGC 147, 
NGC 185, IC 10, LGS 3, and IC 1613 (though not included 
in fitting the plane) within 15° from the pole for both data- 
sets. Within the 30° region, we additionally identify M32 
and And II. For the MI data. And II is also found within 
the 15° region. Both data-sets, MI and KG, thus suggest the 
same satellites to be members of a possible stream, based on 
the intersection of their possible kinematical poles coincid- 
ing with the direction of the normal to the fitted plane. If 
this clustering is true the satellite disc would be rotationally 
supported. The nine satellites are exactly those which Mc- 
Connachie & Irwin (2006b) identified as members of the pos- 
sible candidate streams, namely (iii), (iv), and (v) in their 
paper. In total they identified five further possible streams. 
However, the distance uncertainties for Andromeda and its 
satellite galaxies are so large that an identification of possi- 
ble streams based on the intersection-points alone may not 
be sufficient. With the additional argument of a spatially 
plane-like distributed satellite sample a stronger hint for 
possible streams emerges. 

And VI, which was excluded from the plane fitting as 
an outlier, is the only satellite for which the polar paths can 
be restricted in both data-sets. The resulting RPP is nearly 
perfectly a meridian at Zm3i ~ 170°. Thus, for And VI the 
pole of the angular momentum vector is far off the pole of the 
fitted plane, and the orbit of the satellite can not be located 
within the plane. For the MI data-set the orbit of And V 
can also be restricted and the resulting RPP is very close to 
that of And VI. This may be an indication for a common 
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Figure 9. Polar paths and restricted polar paths (marked with thicker lines) of the Andromeda satellite galaxies for the MI data-set. 
An overall view is shown using an Aitoff-projection. The top-left and bottom-right plots show a Lambert-projection of the northern and 
southern hemisphere, respectively. Solid paths are for the innermost dozen satellites entering the plane-fitting (§ 3.2) while dashed paths 
are the polar paths for those satellites not used. The direction of the angular momentum vector of M33, as derived from measured radial 
velocities and proper motion, is marked by the triangle. Regions of 15° and 30° distance from the fitted pole (Table 2, entry 'ALS') are 
indicated by the closed, thick solid loops. 



direction of the angular momentum vectors of And V and 
VI. 

M33 is the only Andromeda companion for which a 
measured proper motion is available (Brunthaler et al. 
2005). While M33 is located spatially close to the disc of the 
morphological subsample (mss8, §3.3 & Fig. 7), it cannot 
orbit within this disc: if this disc is rotationally supported, 
then the direction of all orbital angular momentum vectors 
must be close to Zmsi = 168°, 6m31 ~ —27°: we restricted 
the polar paths of And VI (and also And V in the MI 
data), members of the mss8 subsample, to this area on the 



Andromeda sky. The direction of the M33 pole is far off 
this position (> 120°), leading to the conclusion that M33 
and the mss8 subsample of satellites cannot have a common 
orbital plane. 

As a caveat we note that the RPP criterion applied 
above to the polar paths depends on two uncertain matters: 
the relative distance uncertainties and the true proper mo- 
tion of Andromeda. Thus we can only treat the restriction 
criterion as a hint for more plausible regions of polar paths. 
Also note that a plane fitting algorithm and the appearance 




of a clustering of intersection points are not independent. 
The intersection point of the polar paths of two satellites 
is the direction of the normal of the plane containing these 
two satellites and the coordinate origin. 



4.3 Kinematic subsample of Andromeda satellites 

Applying ALS to fit a plane for the kinematically motivated 
subsample of nine satellites given above (M32, NGC 205, 
And I, NGC 147, And II, NGC 185, IC 10, LGS 3, and 
IC 1613), but for the time being excluding IC 1613, the pole 
comes out to be located at (/m3i ~ 69.9°, feivisi ~ —35.2°) 
[(«M3i = 73.6°, 6m3i = -35.0°)] (MI[KG]-data-set), with 
axis ratios c/a = 0.12 and b/a = 0.50 [c/a = 0.15 and 
b/a — 0.71], i.e., a highly oblate (disc-like) configuration 
in both data-sets (Table 2). The pole is very close to the 



pole found for the full sample of twelve satellites within 
the approximate virial radius. In Fig. 5 it is clearly vis- 
ible that also some satellites are excluded here (encircled 
by light gray circles) which lie spatially close to the ini- 
tially fitted plane. Including the very distant (and possi- 
bly bound) satellite IC 1613 we find the pole of the fit- 
ted normal to be located at (/m3i = 74.4°, feM3i = —40.5°) 
[(?M3i = 74.7°, feM3i = -40.6°)], with axis ratios c/a = 0.10 
and b/a = 0.35 [c/a = 0.11 and b/a = 0.43]. The latter 
axis ratio, leading to a more prolate configuration, is totally 
dominated by this one very distant satellite (IC 1613). We 
concentrate our further analysis on the sample of eight satel- 
lites without IC 1613, because it is outside the approximate 
virial radius of M31. We refer to this subsample in Tables 2 
- 4 as kss8 (kinematic subsample of eight satellites). 

Applying the AE test to the kss8 subsample the direc- 
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tion of the principal axis is (^msi = 70.5°, 6m3i = —35.2°) 
[(^M3i ~ 74.2°, foMsi = —35.2°)] with a spherical standard 
distance A^ph = 2.4° [A^ph = 1.6°] (KG [MI] data). The 
location of the principal axis is in good agreement with the 
single fits above, and Asph is for both data-sets significantly 
smaller than for the full data-sets, and also an order of mag- 
nitude smaller than for the mss8 subsample (§ 3.3). 

Performing the bootstrap analysis with 5 000 re- 
samplings for the kinematically motivated subsample yields 
shape parameter 7 = 5.9 and strength parameter ( = 4.7 
[7 = 14.9, ^ = 4.2]. For both data-sets the distribution of 
the directions of the fitted planes of the bootstrapped data 
are strongly concentrated and clustered. The derived spheri- 
cal standard distance is Asph = 9.8 [Asph = 11.5], indicating 
that the systematic effects caused by the distance uncertain- 
ties are smaller than the internal scatter as derived by the 
bootstrapping. 

Even though the morphologically motivated subsam- 
ple (§ 3.3) has a smaller rms-height A than the kinemat- 
ically motivated subsample for the KG-data and thus ap- 
pearing as 'thinner' disc, the bootstrap analysis shows that 
the latter one has a more pronounced planar-like feature. 
The strong clustering is found in both data-sets. We have 
therefore uncovered a sample of eight M31 satellites (nine if 
IC 1613 were included) which span a very pronounced disc- 
of-satelhtes that is probably rotationally supported. 



5 STATISTICAL SIGNIFICANCE OF 
DISC-LIKE DISTRIBUTIONS 

In order to study the possible physical nature of the MW 
satellite system the statistical significance of the observed 
anisotropy given a parent distribution needs to be quanti- 
fied. According to the null-hypothesis, the parent distribu- 
tion ought to be a dark-matter sub-halo distribution which 
may be spherical, oblate, or prolate. To evaluate the sig- 
nificance of planar distributed satellite systems we compare 
the bootstrapped samples of the observed distribution with 
bootstrapped data of random samples from the parent dis- 
tribution. For this we first create spherically isotropic distri- 
butions, where the radial linear probability density is pro- 
portional to p(r) oc r~^ , p — 2 psph(f, ■!?, <f)) cx r"**, g = 4) 
consistent with the radial distribution found for the Milky 
Way (Kroupa et al. 2005) and Andromeda (Koch & Grebel 
2006). Random oblate, prolate, or triaxial ellipsoidal distri- 
butions with axis ratios c/a and b/a are then constructed 
by scaling the components of the random spatial position 
vectors while keeping the volume of the ellipsoid invariant. 
As shown in Eq. (1) the formally expected relative height 
of a spherical distribution is dependent on the minimum 
and maximum radii. Therefore the random samples are set- 
up with the minimum and maximum radii as found for the 
Milky Way (see Table 1). For ellipsoidally distributed ran- 
dom samples the initial values are scaled such that the final 
distribution has the expected minimum and maximum radii. 

As for the observed data, each random sample is boot- 
strapped 10 000 times and we calculate the shape parameter 
7 and the strength parameter C, of the resulting distribu- 
tion of fitted normal vectors. Fig. 11 (central panel) shows 
a contour-plot vs. 7 derived for 100 000 random samples 
from an isotropic distribution [a = b = c) each consist- 
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Figure 11. A contour-plot of strength parameter C, versus shape 
parameter 7 for 100 000 random samples each consisting of eleven 
model satellites with an isotropic (central panel), oblate (c/a = 
0.5, b/a = 1.0, top panel), and prolate (c/a = b/a = 0.5, bottom 
panel) parent distribution, each individually bootstrapped 10 000 
times. The shaded contours show the density distribution of the 
derived parameters, dark being high density. The contour lines 
show the enclosed values with significance levels as labelled. The 
star marks the parameters derived for the MW. 

ing of eleven model satellites, bootstrapped 10* times. As 
can be seen from the shaded regions, which show the den- 
sity distribution, the distribution of normal vectors of boot- 
strapped random samples is not expected to be randomly 
distributed in 7, space. They are typically found to be 
girdled or transitional (7 < 1) and marginally concentrated 
(C < 3), while there is also some fraction of clustered (7 > 1) 
distributions. For oblate parent configurations (Fig. 11, up- 
per panel) much more clustered distributions (7 > 1) re- 
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Figure 12. Approximate contour lines for 1, 2, and 5 per cent 
probability that the innermost eleven Milky Way satellites are 
drawn randomly from a parent distribution with initial axis-ratios 
c/a and b/a. The grey shaded region shows the range of axis ratios 
typically found for dark-matter haloes in numerical simulations. 

suit which are typically also more concentrated (C ^ 3). On 
the other hand, for prolate parent configurations (Fig. 11, 
bottom panel) typically a much higher fraction of girdled 
distributions (7 < 1) results, but also with higher concen- 
tration parameter (C 3). Small concentration parameters 
are mostly found for an isotropic distribution. 

To calculate the significance of an observed distribution 
the joint distribution function -D(7, C) is computed and the 
percentile of bootstrapped random samples is derived for 
which both, the shape parameter and the strength param- 
eter, are larger than found for an observed satellite distri- 
bution, e.g. of the Milky Way (Fig. 11, contour lines). For 
each parameter pair of initial values c/a and b/a we create 
100 000 random satellite samples each consisting of eleven 
satellites. Each of these samples is individually analysed us- 
ing the bootstrap method with 10 000 re-samplings. This 
required a large amount of CPU power and we ran the sim- 
ulation on the computer system of the Argelander-Institute 
and the CIP-pool^ of the physics department. The full run 
took about 7 500 CPU- hours running simultaneously on up 
to 30 PCs from 500 MHz class to 3 GHz class CPUs using 
a distributed computing technique. 

5.1 The Milky Way 

The percentile of models found to have bootstrapped dis- 
tributions more concentrated than for the Milky Way are 
listed in Table 6. Fig. 12 shows approximate contour lines 
for 1, 2, and 5 per cent probability that the MW satellites 
(KTB data-set) are drawn randomly from a parent distribu- 
tion with initial axis-ratios c/a and b/a. The values typically 
obtained for Milky Way sized DM haloes (see e.g., Libeskind 
et al. 2005) are shown by the grey shaded region. 

The same analysis for a sample of twelve satellites is re- 
peated with 20 000 random samples and the derived shape 
parameters are compared with those found for the Milky 

^ http://cip.physlk.uni-bonn.d6 



Table 6. Percentile of bootstrapped random samples for which 
the shape parameter and the strength parameter indicate a more 
concentrated distribution of the normals of bootstrapped satel- 
lites than found for the Milky Way satellite system. Different 
setup distributions arc used with axis ratios c/a (along rows) and 
b/a (along columns). The top table gives the percentile for the 
innermost eleven satellites, the bottom table for the innermost 
twelve satellites including the UMa dwarf satellite candidate. Pa- 
rameter combinations for which c/b < b/a, i.e. which arc triax- 
ial and more oblate, are highlighted by a light grey background 
colour. 



b/a 

c/a 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

1.0 0.5 

0.9 0.5 0.5 

0.8 0.6 0.6 0.6 

0.7 0.6 0.7 0.8 1.0 

0.6 0.7 0.8 1.1 1.5 2.0 

0.5 0.7 0.9 1.5 2.4 3.3 4.1 

0.4 0.5 0.9 2.1 4.1 6.3 

0.3 0.2 0.8 2.7 6.9 12.2 

To oJ 

0.9 0.3 0.4 

0.8 0.4 0.3 0.5 

0.7 0.5 0.4 0.6 0.6 

0.6 0.5 0.5 0.8 1.1 1.5 

0.5 0.4 0.7 1.1 2.0 2.7 3.6 

0.4 0.3 0.7 1.6 5.4 3.2 

0.3 0.1 0.6 2.2 6.2 11.6 



Way satellites including the UMa dwarf galaxy. The result- 
ing fractions are listed in the bottom part of Table 6. Includ- 
ing the most recently discovered dSph in Canes Venatici we 
ran only one test with 13 satellites and 20 000 random sam- 
ples for a spherical setup. For this run we find 0.3 per cent of 
the random samples more concentrated than for the Milky 
Way sample. 

The null-hypothesis that the Milky Way satellites are 
drawn randomly from a spherical or mildly triaxial parent 
population can be excluded at very high statistical signif- 
icance (^ 99.5 per cent, confirming the results of Kroupa 
et al. 2005). With increasing triaxiality, the probability in- 
creases for oblate configurations (which are marked with a 
grey background in Table 6). Prolate configurations are ba- 
sically excluded, except for configurations nearly perfectly 
triaxial, e.g. for c/a = 0.5 and b/a — 0.7, where the proba- 
bility may be of order 1 per cent. Including the UMa dwarf 
galaxy increases the significance of this result (reduces the 
propability) . 

5.2 Andromeda 

For Andromeda, using an appropriate setup, the probabil- 
ity that the satellite distribution is drawn randomly from 
a spherically isotropic parent distribution is already 12 per 
cent for the KG-data. This reflects the fact that we find 
the bootstrapped normals of the M31 satellite system to be 
much less clustered than for the Milky Way. The distribu- 
tion of bootstrapped normals for the MI data-set are even 
less clustered and the probability is thus larger. So the hy- 
pothesis that the M31 satellite system is drawn randomly 
from a spherically isotropic parent distribution can not be 
rejected at present, using the available data. 
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5.3 Statistical significance of M31 subsamples 

Performing the analysis as described above for the mor- 
phologically (§3.3) and kinematically (§4.3) motivated sub- 
samples would not yield the correct significance: the null- 
hypothesis for deriving the statistical significance is that all 
satellite galaxies within a certain radius trace a parent dis- 
tribution, i.e. one assumes all satellite galaxies to be of the 
same origin, namely luminous DM sub-structures. But per- 
forming the analysis as above would imply that the eight 
satellites selected by the morphological or kinematical crite- 
rion are luminous DM sub-structures and the excluded ones 
are of another origin. Even in the case of the morpholog- 
ically motivated subsample, where one may argue that all 
dSphs/dEs have been build due to the same mechanism, 
some were excluded because of their large distance from the 
fitted plane. But they should have been included to derive 
the significance, because otherwise one implies a different 
origin of the excluded satellites. Koch & Grebel (2006) did 
not take this into account. 

To derive the significance correctly, the procedure is to 
set up twelve model satellites within the appropriate dis- 
tance range from M31 and then select all combinations of 
eight model galaxies out of the full sample. For each of these 
subsamples a full bootstrap analysis has to be performed. 
Then the fraction of random samples that have subsamples 
of eight model satellites with a more pronounced disc-like 
distribution than the observed sample needs to be calcu- 
lated. 

Practically, this analysis can only be performed for a 
few specific set-ups because the runs take a large amount of 
CPU time. We performed one run for each of the data-sets, 
MI and KG, respectively, using an isotropic parent distri- 
bution. The differences between both random runs are the 
minimum and maximum radii, chosen to match the values of 
the corresponding observed data-sets. We set-up 10 000 ran- 
dom samples. The full run took about 20 000 CPU hours, 
running simultaneously on more than 30 standard PCs. 

The probability that the morphologically motivated 
subsample is picked from samples that are randomly drawn 
from an isotropic distribution is 100 per cent for the MI 
data-set and 95 per cent for the KG data-set. The probabil- 
ities for the kinematically motivated subsample are 17 per 
cent for the MI data-set and 10 per cent for the KG data- 
set. Even though appearing as a 'thin' disc of satellites, we 
conclude that the morphologically motivated subsample is 
just a chance alignment of galaxies. For the kss8 subsample 
we can also not reject the hypothesis that the satellites are 
picked from a random sample, but the probability is much 
lower (< 17%) for both data-sets. 



6 CONCLUDING REMARKS 

We introduce a framework of mathematical methods to anal- 
yse the spatial distribution of the satellite galaxies of the 
Milky Way and Andromeda. A statistical analysis based on 
the 'thickness' of a planar-like distribution alone is not suf- 
ficient to characterise the distribution (Kang et al. 2005; 
Zentner et al. 2005). The bootstrap method is used to derive 
the distribution of poles of fitted planes which are analysed 
using methods based on the statistical analysis of spherical 



data. These methods quantify the robustness of a planar- 
like distribution. Thus, a population of satellites that is not 
planar-like will in our analysis be robustly identified as an 
unstable distribution of poles. We make our implementations 
of the algorithms available under an open source licence at 
http://www.astro.uiii-bonn.de/downloads. Note that the 
analysis presented here is based on more than three CPU 
years on state-of-the-art PCs using a distributed computing 
technique. 

6.1 The Milky Way system 

Applying two methods (ALS and ODR) to fit planes, the 
Milky Way satellite system within 254 kpc is found to be 
highly anisotropical. All companion galaxies are aligned in 
a disc-like structure with a rms-height of only A = 18.5 kpc 
- 22.8 kpc (Table 2). This disc-of-satellites is highly inclined 
with respect to the MW disc, |6mw! ~ 12"^, passing the 
Galactic plane close to the Galactic Centre, Dp ~ 8 kpc < 
A. 

Satellite galaxies are believed to be luminous dark mat- 
ter subhaloes. In recent cosmological simulations MW sized 
DM haloes are found to be triaxial, c/a = 0.6-0.8, typically 
more prolate than oblate (e.g. Jing & Suto 2002; Bullock 
2002; Libeskind et al. 2005; Zentner et al. 2005). Libeskind 
et al. (2005) showed that the distribution of a large num- 
ber of dark-matter subhaloes within these host haloes has a 
similar shape, i.e. the distribution of dark-matter subhaloes 
is a fair tracer of their host halo. This distribution, the large 
sample of subhaloes within the host halo, is the parent dis- 
tribution of the Milky Way satellites if the satellite galaxies 
are luminous, dark matter dominated subhaloes. 

Applying our new statistical framework to the satellite 
system of the Milky Way we show that the hypothesis that 
the MW satellites are drawn randomly from an isotropic or 
mildly triaxial distribution can be excluded at a high sta- 
tistical significance level (^ 99.5 per cent). It can also be 
excluded that the parent distribution has a prolate shape 
as derived from CDM theory. The null-hypothesis that the 
satellite system is drawn randomly from a dark-matter par- 
ent distribution cannot be rejected only if the parent dis- 
tribution is highly triaxial and oblate, i.e. if the parent dis- 
tribution is already disc-like. In this case, and as long as 
no host-galaxy formation is included in large scale CDM 
simulations, the disc of the Milky Way has to be postulated 
to be nearly perpendicularly oriented to the highly oblate 
host halo, because we find the disc-like structure of satellite 
galaxies to have a polar alignment. Even so, the required 
highly triaxial oblate DM-host shape is not consistent with 
the results of modern CDM structure formation simulations. 
Furthermore, recent measurements of the shape of the MW 
potential using the Sagittarius stream show it to be spherical 
within about 60 kpc (Fellhauer et al. 2006). 

The recent discoveries of two additional faint MW com- 
panions increases the confidence of the above statements. 
It therefore follows that any ACDM sub-structure distribu- 
tion is inconsistent with the observed morphology of the 
MW satellite population (Fig. 12). However, since the SDSS 
(York et al. 2000) mostly covers the north pole region of 
the Galactic sky (Fig. 1), newly detected dwarf galaxies are 
very likely close to the fitted disc. To get an answer beyond 
this possible bias it will be crucial to extend the search for 
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MW companions over a larger area of the sky and particu- 
larly at lower galactic latitudes as is planned by the Stromlo 
Missing Satellite Survey (Jerjen et al. 2006, in preparation) 
using the new ANU SkyMapper telescope (Schmidt et al. 
2005). We note though that if additional very faint dwarf 
galaxies are discovered to not lie within the great disc of 
satellite galaxies as quantified here, we are nevertheless left 
with the fact that the eleven most luminous MW satellite 
galaxies are aligned in the disc. 

For the Milky Way a significant fraction of dwarf galax- 
ies may be invisible in the optical due to obscuration by the 
Galactic disc at low galactic latitudes (e.g., Mateo 1998). 
Within the virial radius of the Milky Way about half of 
the total volume has latitude b ^ 30°. Andromeda is in 
that sense a better probe since it's halo is not that much 
affected by obscuration (McConnachie & Irwin 2006b). A 
simple estimate for the MW, assuming that all undetected 
satellites with b ^ 15° are obscured, 50 per cent of all un- 
detected satellites with 15° < b ^ 30° are obscured, and 
assuming that the undetected satellites are homogeneously 
distributed over the whole sky, we find that about 35 per 
cent of all satellites with b ^ 30° may be found more than 
la — A and about 30 per cent more than Scr off the fitted 
disc-of-satellites . 

6.2 The Andromeda system 

For the Andromeda satellite system it cannot be excluded 
that it has been drawn randomly from a spherical isotropic 
parent distribution. However, we do find the M31 satellite 
system to be anisotropic, but the details depend on the data- 
set used. The fitted disc-like structure for all satellites within 
the approximate virial radius is not as polar-aligned as for 
the MW (|foM3i| ~ 30°) and it is approximately twice as 
'thick' as found for the MW (Table 2). 

Two incompatible subsamples of satellite galaxies which 
have a disc-like distribution can be identified: one morpho- 
logically motivated as proposed by Koch & Grebel (2006) 
and one kinematically motivated (§4.3). Since the disc-like 
satellite system of the MW is dominated by dSph galaxies, 
one can speculate about a common building mechanism for 
all the dSphs in the Local Group. If this mechanism was 
the break-up of a large, gas-rich galaxy or the formation of 
tidal-dwarf galaxies (TDGs) in an early major-merger event, 
one expects the dSphs to have initially correlated directions 
of their angular momentum vectors, supporting a disc-like 
structure. This would favour the morphologically motivated 
subsample since it is initially build up of dSphs only. How- 
ever, at least two of the Andromeda dwarf spheroidals within 
the virial radius do not fit into this picture and have to be 
excluded because of their apparently large distance from the 
plane. Also only one out of three morphologically similar dEs 
is found close to the disc. But we cannot exclude the possi- 
bility that some massive TDGs may retain their interstellar 
media to appear today as dirr galaxies (Hunter et al. 2000; 
Recchi et al. 2006). 

The disc-like distribution of the dSph/dE galaxies 
around Andromeda is present for the KG data-set only, but 
can not be identified for the MI data-set. Comparing the 
results of the AE and the bootstrap test it follows that even 
in the KG data-set the systematic errors caused by the dis- 
tance uncertainties are larger than the intrinsic scatter as 



quantified by the bootstrapping. The high significance of the 
morphologically motivated disc as derived by Koch & Grebel 
(2006) is incorrect because of three reasons: firstly they de- 
rived their significance based on the thickness of the disc 
alone, secondly they used the ODR fitting method which is 
affected by the systematic alignments of the distance uncer- 
tainties of the M31 satellites, and thirdly they did not use 
the correct null-hypothesis to derive the significance. We 
find the distribution of the morphologically motivated sub- 
sample of M31 satellite galaxies to be fully consistent with 
being picked from a random distribution. 

Combining the polar-paths of the Andromeda satel- 
lites with the fitted plane it follows that there may be a 
kinematic association of some of the M31 satellites: M32, 
NGC 205, And I, NGC 147, And II, NGC 185, IC 10, LGS 3, 
and probably, though due to its large distance unlikely, 
IC 1613. These satellite galaxies were also identified as possi- 
ble stream members by McConnachie & Irwin (2006b), while 
there are also other possible streams which now seem to be 
less likely. This sample of eight Andromeda satellites forms 
a very pronounced thin disc with inclination of ~ 31° away 
from a polar alignment, and that holds true for both data- 
sets (MI and KG). The kinematically motivated subsample 
is found to have a much more pronounced disc-like distribu- 
tion than the morphologically motivated one, although its 
thickness is larger. Even though the kinematically motivated 
subsample has a much higher statistical significance than the 
morphologically motivated one, we can also not exclude the 
possibility that this sample is picked from a random distri- 
bution. 

We argue that the close proximity of M33 to the disc of 
the morphologically motivated subsample is a pure chance 
alignment since the direction of its angular momentum vec- 
tor can not be correlated with that of the M31 dSphs. Using 
the RPPs (§4.1) we show that, if the dSphs have a common 
direction of their angular momentum vectors, it is far ofi^ 
(> 120°) the pole of the M33 orbit as derived from its mea- 
sured proper motion. But M33's kinematical pole lies close 
to the pole of the whole distribution and to the pole of the 
kinematically motivated subsample (Figs. 9 and 10). 

Using radial velocity measurements we created re- 
stricted polar paths, showing that at least the direction of 
the angular momentum vector of And VI can be restricted to 
an arc segment of 180° , i.e. the sense of rotation of the orbit 
around Andromeda can be restricted in both data-sets. For 
six of the M31 satellites the possible poles can only be re- 
stricted for one of the data-sets. Restricted polar paths hint 
to more likely regions of the locations of angular momentum 
vectors on the M31 sky. 

The global infiuence of M32 and M33 on the satellite 
system of Andromeda remains to be studied. The Triangu- 
lum Galaxy M33 may be sufficiently massive (> 5x W^^ Mq, 
Herrmann & CiarduUo 2005) to mix up the satellite system 
of Andromeda or bring its own satellite galaxies into the 
M31 system. Similarly the cE M32 may have had a strong 
infiuence on the dynamical evolution of Andromeda's satel- 
lite galaxies if massive enough in the past (e.g., Bekki et al. 
2001). 
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Figure 13. An AitofT projection in supergalactic coordinates of 
the directions of the spin poles of the MW and M31 discs marked 
by filled circles. The directions of the normals of the fitted plane 
of all satellites of the MW and M31 within their virial radii are 
marked by squares and triangles, respectively (MI data-set in dark 
grey, KG data-set in light grey). 



6.3 The Milky Way vs. Andromeda 

There is no obvious evidence for a spatial association of the 
disc-like structure of the Milky Way satellites and that of 
the Andromeda satellites. If M31 (located at Zmw = 121.7°, 
6mw = —21.5°) and its satellites were associated with the 
disc-of-satellites of the MW, M31 ought to be close to the 
fitted disc-like structure and the fitted planes ought to be 
aligned. Instead, M31 is « 55° off the fitted plane of the 
MW satellites. Similarly, the angle between the normals of 
the fitted satellite planes of the MW and M31 is « 50° - 
60°, but would be 0° if they were perfectly aligned. 

Fig. 13 shows the directions of the spin-poles of the 
MW and M31, as well as the directions of the normals of 
the fitted disc-like structures of their satellite distributions 
in supergalactic coordinates (de Vaucouleurs et al. 1991). 
The directions of the spin-poles (SP) are marked by filled 
circles. The directions of the normal of the fitted plane to 
the MW satellite system is marked by squares, the filled 
one pointing in the direction that is close to the kinematical 
poles as derived for some of the MW satellites (Palma et al. 
2002, and references therein; Kroupa et al. 2005) indicating 
the sense of rotation if the disc-like structure is rotation- 
ally supported. The direction of the fitted normal to the 
M31 satellite system marked by triangles is arbitrary ('up' 
or 'down'). It follows that both spin-axes of the disc galax- 
ies lie within about 30° of the supergalactic plane and that 
the disc-of-satellites of both hosts are likewise oriented such 
that their poles lie within 30° of the supergalactic plane. All 
discs are thus highly inclined with respect to the supergalac- 
tic plane, but do not appear to be mutually aligned. This 
would appear to contradict the notion that the satellites ac- 
creted individually preferentially from the direction of the 
supergalactic plane, i.e. from the direction of the medium- 
scale (~ few Mpc) matter distribution. 

The highly-inclined orientation of the stellar discs of 
the MW and of M31 relative to the supergalactic plane 
can be understood either as being a result of tidal torquing 
(Navarro et al. 2004) or resulting from the perpendicular col- 
lapse of matter onto the supergalactic plane (Doroshkevich 
1973; Doroshkevich et al. 1978; see e.g. also Hu et al. 1998). 
The disc-of-satellites of the MW and of M31 are both also 
highly inclined to the supergalactic plane and the respec- 



tive galactic discs of their hosts. Tidal torquing, if respon- 
sible for the orientation of the galactic discs, can therefore 
not be the origin of the disc-of-satellite orientations. Instead 
these two discs could also result from the collapse of matter 
onto the supergalactic plane. However, this would beg the 
question as to why they are so highly inclined to the host 
galactic discs. Alternatively, it would appear more natural 
or intuitive to understand the disc-of-satellites as being the 
result of stochastically occurring mergers which leave pop- 
ulations of related TDGs. Such populations would remain 
visible for highly inclined events relative to the host discs, 
because populations of TDGs in low-inclination orbits would 
precess apart and possibly end up merging altogether with 
their host discs (Penarrubia et al. 2002). 
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APPENDIX A: CONDITIONAL 
BOOTSTRAPPING 

To uniquely perform a plane fitting, at least three different 
objects in a bootstrap sample are required. To correctly cal- 
culate the total number of possible distinct bootstrap sam- 
ples, this can generally be formulated as: 

What is the number of ways of picking k unordered out- 
comes with replacement from n possibilities under the con- 
dition that at least q outcomes are different? 

Without the condition that at least q outcomes are dif- 
ferent, the total number is given by the number of unordered 
samples with replacement 



iVtot = 



k-1 
k 



(Al) 



We start by calculating the number of ways of pick- 
ing exactly p — q different outcomes. This is given by the 
number of ways of picking p unordered outcomes without 
replacement from n possibilities times the number of ways 
of picking k — p unordered outcomes with replacement from 
p possibilities. 



p + (fc — p) — 1 
k — p 



k-1^ 
k-pj 



(A2) 



This can be understood as follows: 



1 I ^'2 



e {ji,j2,...,jp}. 



First choose exactly p different outcomes from all ?i pos- 
sibilities without replacement (all p must be 
different); the number of ways is given by (^) . Since exactly 
p different outcomes were premised, it can now only be cho- 
sen from those p outcomes that have been chosen in the 
first step (so now the number of possibilities is p) , which are 
r'lZ''-') ways. 
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Next calculate this for p 
add up all numbers of ways: 



1 and so on and finally 



E 



k-i 

k — p 



(A3) 



Since the total number of possible distinct bootstrap sam- 
ples is typically very large (e.g. ^A'tot = 77 557 275 for 
n = 15) a sufficiently large number of re-samplings is used, 
as is usually done in bootstrapping. 

Note that bootstrapping is different to the sampling 
method Koch & Grebel (2006) employed, while mislead- 
ingly referring to it as bootstrap. These authors fitted a 
plane through every possible combination of 3, 4, . . . , n ob- 
jects, an approach that is related to the jackknife method. 
The different sampling method may have biased the results 
of Koch & Grebel (2006). 



APPENDIX B: EXTRASOLAR CELESTIAL 
COORDINATE SYSTEMS 

In order to study the three-dimensional distribution of a 
satellite system it is most convenient to transform their po- 
sition vectors relative to the observer into an extrasolar co- 
ordinate system. We elaborate the transformation in much 
detail in the most general way, using vector and matrix op- 
erations, to allow for an as wide as possible range of applica- 
tions. The algorithm described also allows for a free choice 
of the reference point of the / = 0° direction of the target 
system. Similar transformations (in a less general way) were 
used in two recent relevant works on the spatial distribution 
of the Andromeda satellite galaxy system (McConnachie & 
Irwin 2006b; Koch & Grebel 2006), but the transformations 
differ significantly. As we will show the transformation as 
applied by Koch & Grebel (2006) is incorrect. 

The basis system for the transformation may be any 
celestial coordinate system, e.g., equatorial coordinates (a, 
S) or Galactic coordinates {I, b). The origin and the orien- 
tation of the target coordinate system must be given as (i) 
positional coordinates, (ii) inclination, and (iii) position an- 
gle in the basis system. We exemplify the transformation by 
transforming into an Andromeda-centric coordinate system. 
In the following, bold symbols like r denote vectors, whereas 
r is the absolute value of a vector r = |r|. Unit vectors in a 
specific direction are denoted as — r/r. 



Bl An Andromeda-centric coordinate system 



Let (e^, , fiy , ) be a Cartesian coordinate system cen- 
tred on Andromeda similar to the galactocentric coordinate 
system (el^", e^™, e^™) of the Milky Way (see also Mc- 
Connachie & Irwin 2006b). We define the coordinate system 
such that e"^^ is the direction of the normal to the stellar 
disc of Andromeda, e^'^^ is the projected direction from the 
MW to M31 onto the disc of M31 and e^^' = ef^^ x e"^\ 
Equatorial coordinates can be calculated as follows: 



, y 

arctan | — 



— arcsm 



(;) 



(Bl) 
(B2) 
(B3) 




Observer 



Figure Bl. This artist- view shows the orientation of the normal 
n to the equatorial plane of Andromeda which is used to calculate 
positions in the Andromeda-centric coordinate system, is the 
direction from the observer to Andromeda. Cp is parallel to the 
celestial equator and pointing towards east and is pointing 
towards the celestial north pole. The inclination i is the angle 
between the normal of the equatorial plane and the line-of-sight. 
The position angle 6 is the angle between the projected major axis 
of the disc of Andromeda and the direction to the north celestial 
pole, measured from north through east. 



where x,y, z are the components of a vector either given 
in the galactocentric or in the Andromeda-centric Cartesian 
coordinate system. Note that if we give galactocentric coor- 
dinates (/mw,6mw), the centre of the coordinate system is 
the Galactic Centre (GC), whereas the centre of the stan- 
dard 'Galactic coordinate system' {l,b) is the Sun. 

To calculate the Cartesian components of a vector in 
the Andromeda-centric coordinate system, the heliocentric 
position vector r is shifted to the centre of M31 and rotated 
to the Andromeda-centric coordinate system. 



r-Msi = Rm3i {r - r(M3l)) , 



(B4) 



where r(M3l) is the heliocentric position vector of M31. 

To calculate the rotation matrix Rmsi in Eq. (B4) it is 
not necessary and also not advisable to start with coordi- 
nates given in Galactic coordinates {l,b). Rather, it is rec- 
ommended to work from right ascension and declination (a, 
S). The reason is the way the orientation of the Andromeda 
disc is given: the position angle 9 is the angle between a 
given direction and the direction to the north celestial pole 
(NCP), conventionally measured from north through east 
(e.g., Binney & Merryfield 1998). Using equatorial coordi- 
nates eliminates the need to correct for the difference be- 
tween the direction to the NCP and the north Galactic pole 
(NGP) which is in general different for different positions on 
the celestial sky. 

The first step is to calculate the transformation-matrix 
that rotates the 'normal triad' (Murray 1983) parallel to the 
axes of the Cartesian coordinate system. This is given by the 
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Br, Cp, and Bq 

Rrpq{o:,S) 



(fir 



Bp Bq) 

''cos S COS a 
cos 5 sin a 
>. sin 5 



— sin a — sin & cos a\ 
cos a — sin (5 sin Q 1 (B5) 
cos 5 / 



Br is in general pointing in the direction with equatorial co- 
ordinates (a, 5), fip is parallel to the celestial equator, posi- 
tive towards the east, and e, is pointing towards the north 
(Murray 1983). Here (e^ Bq) is the normal triad of the 
coordinates of M31 (qmsi, <5m3i)- 

A preliminary rotation matrix is given by: 

Rm31 = Ry(90° - i)Ra:(90° - R.pq (qm3i , (5m3i ) (B6) 

where R^^yj ^^"i) are the matrices which perform a counter- 
clockwise rotation about the Cartesian coordinate axes by 
an angle 7, e.g. 





COS 7 

— sin 7 
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sin 7 
cos 7 1 



(B7) 



Ri(90° — Q) corrects for the position angle and Ry(90° — i) 
for the inclination of M31 (Fig. Bl). This rotation is pre- 
liminary because it performs a transformation such that the 
Sun is located at (/M3i,feM3i) = (180°, -12.5°) but we de- 
fined the Milky Way, i.e. the Galactic Centre, to be located 
at Im31 ~ 0°. Therefore an additional rotation of the order 
arctan(8.5 kpc/785 kpc) = 0.6° hcis to be performed. 

The procedure to calculate this angle is as follows: 
(i) transform the vector pointing from Andromeda to the 
GC to the preliminary system 



r'(GC) = Rm31 (r(M3l) - r(GC)) 
(ii) project this vector onto the plane of M31, 

r'(GC))-fe'r( 



/M31 / /M31 

e y ■ 



(B8) 



r'(GC)), 
(B9) 



and (iii) calculate the angle /3 between this vector and the 



. /M31 

prelimmary e ^ -axis: 
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Now the full transformation matrix can be calculated: 



R?(/3)Rj,(90° -i)R490° -e) 

RT-pq(aM31, <5m3i) 



^ 0.7703 
-0.6321 
,-0.0839 



0.3244 
0.5017 
-0.8019 



0.5490N 

0.5905 

0.5915; 



(B12) 
(B13) 



Note that Rj^(/3) = R,(— /3). We use an inclination i = 
77.5° and a position angle = 37.7° (de Vaucouleurs 1958). 
To calculate the angle j3 and the matrix Rrpq{aM3i, Smsi) 
in Eq. (B13) we use a distance of 785 kpc and coordinates 
aM3i = 00''42'"44.3, Smai = -h41°16'09.4" for M31, and we 
employ a distance 8.5 kpc of the Sun from the GC. 

The coordinate transformation as given in Koch & 



Grebel (2006) must be wrong as the orientation of the triad 
at the position of Andromeda was ignored, just shifting to 
the centre of M31. The offset of the directions of the NCP 
and the NGP was also not correctly calculated. This offset 
is by chance very small at the position of Andromeda on 
the celestial sky (compare to McConnachie & Irwin 2006b). 
According to their transformation, Andromeda would have 
an inclination of i = 66° and a position angle of 6^ = —16°. 
However, since they performed only orthogonal transforma- 
tions, translation and rotation, the relative orientations of 
the satellites are preserved. 



B2 General transformation 

Equation (B4) is the transformation equation performing 
the necessary linear operations, translation and rotation, re- 
spectively. Generally this becomes 



Ri 



(rs - r8(Or)) , 



(B14) 



where rs is the Cartesian vector of any object picked to be 
transformed to the target system T, rs(OT) is the vector 
of the origin of the target system T, both vectors given in 
the basis coordinate system B. The appropriate rotation ma- 
trix Re_>T can be calculated as given in Eq. (B6) or (B12). 
Equation (B6) gives the rotation matrix for which the Sun is 
located at longitude 180°. Equations (B8) - (B12) show how 
to derive the rotation matrix such that the GC is at longi- 
tude 180°. rr is the Cartesian vector of the object given in 
the target coordinate system. Equatorial coordinates follow 
from Eqs. (Bl) - (B3). 



B3 Example Transformation 

We give here an example of how to transform the coordi- 
nates of a companion galaxy into the Andromeda-centric 
coordinate system as described in Section Bl. We apply 
the transformation to M33 with coordinates a = 0l''33™51, 
S = 30°39'36", and a distance D = 809 kpc from the Sun, for 
Andromeda, a = 00''42™44.3, S = +41° 16'09.4", with a dis- 
tance D — 785 kpc. First we need to calculate the Cartesian 
coordinates of M31 and M33 in the heliocentric coordinate 
system (Murray 1983): 

r(M3l) = (579.791 109.391 517.785) , 
r(M33) = (638.372 277.075 412.543) , 

where ^ denotes transposition. Next we apply Eq. (B4) using 
the matrix Rmsi as given in Eq. (B13): 

rM3i = RM3i(r-(M33) - r(M3l)) , 

= (41.737 -15.051 -201.635) , 

and finally we can calculate the galactic coordinates in the 
Andromeda-centric coordinate system (Eqs. B1-B3): 

/m31 = 340.2° , 
6m3i = —77.6° , 
'"M31 = 206.5 kpc . 
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APPENDIX C: SOFTWARE 
IMPLEMENTATION 

All algorithms are implemented using the Python program- 
ming language (http : / /www . python . org) and extending it 
using C routines. We use the packages Numeric (http: 
/ /numpy . scipy . org) and its implementations of basic linear 
algebra routines as well as the SciPy package (http : / /www . 
scipy.org). The ODR method (http://www.netlib.org/ 
odrpack) is implemented using the wrapper-package odr by 
Robert Kern (http://www.python.net/crew/kernr/). 

The software is available at http : //www. astro . 
imi-boim.de/downloads released under the GNU General 
Public Licence (GPL). 
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